Tensor power spectrum with holonomy corrections in LQC 
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In this paper we consider tensor perturbations produced at a bounce phase in presence of the 
holonomy corrections. Here bounce phase and holonomy corrections originate from Loop Quantum 
Cosmology. We re-derive formulas for the corrections for the model with a scalar field content. 
Background dynamics with a free scalar field and multi-fluid potential are considered. Since the 
considerations are semi-classical effects of quantum fluctuations of the background dynamics are not 
taken into account. Quantum and classical back-reaction effects are also neglected. To find spectrum 
of the gravitational waves both analytical approximations as well as numerical investigations were 
performed. We have found analytical solutions on super-horizontal and sub-horizontal regimes 
and derived corresponding tensor power spectra. Also occupation number rik and parameter figw 
were derived in sub-horizontal limit, leading to its extremely low present value. Final results are 
numerical power spectra of the gravitational waves produced in presence of quantum holonomy 
corrections. The obtained spectrum has two UV and IR branches where Vt oc , however with 
the different prefactors. Spectrum connecting these regions is in the form of oscillations. We have 
found good agreement between numerical spectrum and this obtained from the analytical model. 
Obtained spectrum can be directly applied as initial conditions for the inflationary modes. We 
mention possible resulting observational features of the CMB in particular B-type polarization. 



I. INTRODUCTION 



In the Minkowski background free gravitational waves 
fulfil the wave equation {d^ — V^)ft.^^ = 0. Solutions of 
this equation are plane waves /i^i/ = P^^e*''^ ''^''''*', here 
P^i, is polarization tensor. However when cosmological 
expansion is turned on (we assume flat FRW background 
here) additional term appears and the equation of mo- 
tion is modified to {d^ + 3Hdt — V^)^-^^ = 0, where H is 
Hubble parameter. We see that cosmological term acts 
as effective friction. When Universe undergoes expansion 
then H > and gravitational waves are damped. This 
situation corresponds to the present stage of the evo- 
lution. However when Universe is in contracting phase 
H < friction term becomes negative and gravitational 
waves are amplified. Such a phase of contraction is gen- 
eral prediction of Loop Quantum Cosmology (LQC) 3. 
Contracting phase appears also in the string theory based 
theories of the Universe 0, Q and many others. However 
in the present paper we concentrate on the LQC mod- 
els where contracting and expanding regimes are joined 
by the phase of bounce [4]. In the last years perturba- 
tions during the bounce phase were studied extensively. 
Recent review on this issue can be found in ^] . In partic- 
ular perturbations in the Quintom Bounce were studied 

in ig. 

In LQC many physical results can be traced in semi- 
classical approximation. In particular dynamics of the 
Universe can be recovered from the quantum corrected 
Friedmann equation [1, . Similar approach can be also 
applied to describe quantum gravity effects on pertur- 
bations [1ft] . in particular to gravitational waves. It is 
however worth to stress that such an approach is rather 



'Electronic address: ' jakub. mielczarek@uj.edu.pl| 



heuristic and results obtained have to be verified by the 
purely quantum considerations. In particular it has not 
been proved yet whether phase of bounce is generally re- 
alized for the inhomogeneous loop cosmologies. However 
some recent studies show that in case of loop quantized 
inhomogeneous Gowdy spacetime, singularity is avoided 
[nj . In our approach inhomogeneities are treated pertur- 
batively and we neglect their back-reaction on the back- 
ground dynamics. However in the more detailed studies 
these effects have to be also taken into account. In the 
semi-classical approach applied here, quantum gravity ef- 
fects are introduced by the corrections to the classical 
equations of motion. For tensor modes in LQC these 
effects were preliminary studied in Ref. [13, [3] ■ Later 
improved approach was developed introducing holon- 
omy corrections. Results of this paper are a backbone of 
our investigations. In this paper we assume that these 
corrections are valid during the whole evolution. Some 
preliminary studies of influence of the holonomy correc- 
tions for the gravitational waves production have been 
done [l^ [l^, |l3 • However in that papers effects of the 
corrections to the source term were neglected. While 
in classical approach this term vanishes (within linear 
regime) in quantum regime it does contribute. In the 
present paper we improve these studies including a source 
term. 

Besides the holonomy corrections also inverse volume 
corrections are predicted in the framework of LQC. Ef- 
fects of inverse volume corrections on gravitational waves 
were recently studied in Ref. [H, [l^. However in the 
flat FRW background inverse volume corrections exhibit 
flducial cell dependence. This makes those effects harder 
to interpret. However, in the curved backgrounds this 
problem disappears. Since holonomy and inverse volume 
effects differ qualitatively, they should be studied sepa- 
rately. In this paper we follow this line of reasoning. We 
consider consistent model where holonomy corrections in- 
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fluence both background and perturbations parts. 

The organisation of the text is the foUowing. In sec- 
tion |TT] we introduce the equation for tensor modes with 
holonomy corrections. Then in section|TTT]we define back- 
ground dynamics. We consider both the model with free 
scalar field and with multi-fluid potential. Subsequently 
in section IIVI and section |V] we investigate analytically 
and numerically the evolution of the tensor modes. Ef- 
fects of holonomy corrections are investigated. With use 
of numerical computations we calculate power spectra of 
the gravitational waves. In section IVTl we summarize the 
results. Finally in Appendix we introduce gravitational 
waves in LQC framework, derive particular form of the 
holonomy corrections and explain the employed notation. 



II. GRAVITATIONAL WAVES WITH 
HOLONOMY CORRECTIONS 

Equation for tensor modes with LQC holonomy cor- 
rections derived in [14] is given by 



+ 2~k^K - ^^K + TqK = leTrcn^Q., 



(1) 



the quantum holonomy corrections simplify to 



Tc 



n 



Qa 



3 Pc 



1-^ {2V-p)hl 
2 Pc 



(7) 

(8) 



These expressions were first derived in Ref. [T^. How- 
ever we have found a discrepancy between the expression 
for Hq^ derived here and this found in Ref. To 
approve the result presented here we show intermediate 
steps of derivation in Appendix. The difference is 1/3 
factor inside the bracket. To derive these corrections we 
have applied so called /2 scheme of quantisation. Namely 
we used p, = yj /S./p where A = 2\/?,iT^lpy It is well mo- 
tivated to use this particular form of the function [20j . 
However, other choices are in principle also permitted. 
In this paper we consider only /2 scheme, which seems to 
be the best motivated. 

Now equation for the tensor modes ([T]) simplifies to 

^hi + 2k-^K-v^K + fQK^o (9) 

ai]'' drj 
where we have defined the total holonomy correction 



where 



„ „ , p dp\ _2 2 ( smpjk 
Tq = -— Ut 7 — 



1 dH^ 

3Vb dp 



— ^ cos 2^7A; 



SH„ 



(2) 
(3) 



For details and explanation of the employed notation we 
send to Appendix. To derive specific form of the func- 
tions Tq and Hg^, matter content must be defined. In 
this paper we consider models with a scalar field. We 
consider both free and self-interacting fields. In that case 
matter Hamiltonian is up to the second order 



Hm. — 



4 7e Vp V 2 p3 



where homogeneous part is given by 



IK. 



(4) 



(5) 



Here integration was constrained to fiducial volume Vq. 
Further physical results do not depend of this quantity. 
Energy density can be now defined as 



1 dH^ 
Vopy^ dN 



(6) 



Tq^Tq- WttGHq = WnGpy - . (10) 

Therefore also source term correction has been included. 
This is in contrast with the analysis performed in [isl . 
[igI , [TtI . where this influence was neglected. In the clas- 
sical theory in fact this term vanish in the linear order. 
Therefore when fluctuations of vacuum are considered, 
higher order term can be set to zero. However since, due 
to quantum corrections, source term contribute linearly, 
there is no reason to neglect this term. Therefore in the 
present paper we take it into account. 
We introduce new common variable 



aha 



(11) 



where h\ = —h^ = h(Q, h\ = h\ = hf^ and a = y/p. Then 
performing the Fourier transform 



u(?7,x) 



(27r)3 



u(?7, k)e 



ik-x 



(12) 



one can rewrite the equation Q in the form 



(77,k) + [fc^+m^ff]u(r;,k) = 0, (13) 



where /c^ = k • k and 



(14) 



When matter content is defined one can derive partic- 
ular form of the functions ^ and ([3]). Expressions for 



In this paper we aim to solve equation psp . However 
first we must specify the background dynamics. 
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III. BACKGROUND DYNAMICS 



B. Scalar field with a multi-fluid potential 



Background dynamics is governed by the effective 
Friedmann equation 



where 



2pdt) 3^^ 



V3 



(15) 



(16) 



is critical energy density. This equation can be derived 
combining the Hamihon equation p — {p, + iJ^,^''"} 
with the scalar constraint 5„ + H^^"^ = 0. 

Evolution of the scalar field component is governed by 
the Hamilton equations 



dV 



(17) 
(18) 



Energy density and pressure of the homogeneous scalar 
field are expressed as follows 



P4> = I'P^ + Vi^b), 



Free scalar field 



(19) 
(20) 



Energy density of the free scalar field has the form 



1-1 
2 p3 



(21) 



and the solution of the effective Friedmann equation (ef- 
fective background equation) is the following 



where 



(22) 



(23) 



Solution (|22p represents non-singular bouncing evolution 
and is discussed in Ref. 

For further applications, it will be useful to relate the 
coordinate time with the conformal one drj ~ dt/a{t). 
Assuming that rj{t = 0) = 0, we obtain 



riit) 



AV6 



iFi 



113 B 
2' 6' 2'"Z^ 



(24) 



One can shown [21] that the restriction ~ wp^p, 
where w — const in the framework of effective LQC, leads 
to the potential in the form 



cosh 



^67rG(l -l-u;) 



(25) 



Solution of the equations of motion with this potential 
has simple analytic form 

p{t)^Pc{l + 6TTGpc{l + wft^)''^. (26) 
It is worth to mention that for w = 1 and taking 

1 



Pi = A = -K^;7 A 

we recover the solution (|22p . 

In analogy with the free field case we derive 



(27) 



2F1 



1 



1 



2' 3(l + w)' 2' 



(28) 



IV. ANALYTICAL CONSIDERATIONS 

The theory of cosmological creation of particles bases 
on idea of "freezing" of the vacuum fluctuations. On the 
mathematical level this process can be seen as a squeez- 
ing and displacement of the vacuum state |0). This is 
equivalent with the creation of particles. For the non- 
interacting field theories the wave function is a product 
of the functions for the particular modes. Therefore the 
degree of squeezing and coherence can be different for the 
particular modes and is determined by the cosmological 
evolution. The typical scale for which squeezing and dis- 
placement of the vacuum becomes important is the Hub- 
ble scale. Modes of quantum fluctuations becomes classi- 
cal (are described be the coherent states) when crossing 
the Hubble radius. 

To describe process of particles creation quantitatively 
one can consider Bogolyubov transformation between ini- 
tial and final states. Then computing the so-called Bo- 
golyubov coefficients the number of produced particles 
can be obtained. However on the super-horizontal scales 
one can in principle obtain ujI = k"^ + m^fj < and the 
interpretation in terms of particles fails. Then quantum 
state cannot be interpreted in terms of particles. There- 
fore, and from other reasons, it is useful to consider a 
correlation function which is well defined quantity for all 
energy scales. The correlation function for the tensor 
modes takes the form 



|/i?(x,ry)/i^(y,,7)|0) 



f dk^ , sin kr 



kr 



(29) 
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where we have defined the power spectrum 



(30) 



The power spectrum can be related later to the ampli- 
tude of the CMB fluctuations. Therefore it is crucial to 
determinate this function. 

Another way to describe physical properties of the 
quantum state is the mentioned method of Bogolyubov 
coefficients. The relation between annihilation and cre- 
ation operators for the initial and for the final state is 
given by the Bogolyubov transformation 



&k = B+{k)cn, + B^{k)*a\ , 
S+(A:)*4 + S_(fc)a„k , 



(31) 
(32) 



where — — 1- Since we are working in 

the Heisenberg description the vacuum state does not 
change during the evolution. It results that 6k|0in) = 
(fc)*alj.|Oin) is different from zero when B-(k)* is a 
nonzero function. This means that in the final state the 
graviton field considered is no more in the vacuum state 
without particles. The number of produced particles in 
the final state is given by 



|0i„) = |S_(fc)|^ (33) 



The energy density of gravitons is given by 

\B-{k)\ 



Anaj'^dLu /,m2 
dpgw — 2 ■ fiuj ■ \ u 1 1„\\ 



(27rc)3 



(34) 



where we used definition (|33p . To describe the spectrum 
of gravitons it is common to use the parameter 



rigw(i^) = 



P* 



dv 



(35) 



where is the energy density of gravitational waves 
and p.^, is the present critical energy density. 



Free scalar field 



Based on solution ((22l) we derive 



4 {A + Bt^) 



5/3 



> 0. 



(36) 



We show this function in Fig. [TJ We compare it with the 
classical expression 



4 {A + Bt^) 



The difference is significant since now effective mass is 
a non-negative function, rn^g > 0. One can also com- 
pare this with the case when the source term corrections 




FIG. 1: Evolution of the effective masses m^g and rn^^iTQ 
0). In this figure we have assumed = 0.1 Ipi 



were neglected. Then, as can be found in Ref. [l6|, the 
effective mass is negative in some regime and behaves 
like the classical one. Here difference is crucial and has 
important consequences. Namely since TOgg > we have 
cj^ > and interpretation in terms of particles is possible 
on all scales. This becomes a nice property also when the 
Hamiltonian of the perturbations is minimized to find a 
proper vacuum state. It can be shown that when ujj, > 0, 
the Hamiltonian has a minimum for all k and a well de- 
fined vacuum can be found. Otherwise for some k < fcx, 
the lowest-energy instantaneous vacuum state does not 
exist. 

Now we are going to consider the pre-bounce limit. 
Taking \t\ — > oo, we find 



'off 



1 1 
4 77^' 



(38) 



The normalised solution of the equation has the form 



^/4 4fV^^o"(-^fc)- (39) 
V 2k 



u{k,ri) 



We have chosen here advanced modes and performed nor- 
malisation with use of the Wronskian condition. In the 
super-horizontal limit —rjk ^ 1 we can apply the approx- 
imation 



1 + .^ 



In 



(f) 



IE 



(40) 



where — 0.57721 ... is Euler-Mascheroni con- 

stant. Expression for the power spectrum in the super- 
horizontal limit is therefore 



Vrik) = Ak^ { 1 



In - 



ki] 



where 



A 



2 ^3A'/' flpi 



(41) 



(42) 
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To investigate k dependence in formula (|4ip we define 
the spectral index 



riT = 



dlnVTik) 
dink 



and obtain 



In - 



- IE 



1 + ^ 



ln(-^? 



IE 



(43) 



(44) 



We show this function for some fixed time in Fig. [21 
We find that the resulting spectral index is blue and ap- 




FIG. 2: Running spectral index on the super-horizontal scales. 

proaching nx = 3 for fc — > 0. This blue-tilted spectrum 
was predicted earlier in [16]. Recent investigations sug- 
gest that also for inflationary cosmology with holonomy 
corrections obtained spectrum is blue-tilted and tit — 3 
at super-horizontal scales [l7j . 



B. Multi- fluid potential 



where 



9 (l-w)^ 



(49) 



4 (1 + 3w)2' 

Power spectrum of the perturbations is then given as 

7'T(fc) oc (-fc?7)3-2IH^ (50) 

where super-horizontal approximation 
i7W(x)^ 



-r(n) 



(I)"" 



for X ^ 1 



(51) 



has been used. It must be stressed that the above formula 
does not hold for n — (w=l). In that case another 
expansion (|40p must be applied. 

It is worth to mention that scale invariant spectrum 
= I is recovered both for w — —1 and w = 0, as 
it can be directly seen from ([^5]) . This duality was in- 
vestigated in Ref. [l^] in context of the free scalar field 
perturbations. 



C. Sub-horizontal solutions 

Since now we were only concerned with the prc-bounce 
phase. Now we are going to evolve modes through the 
bounce. We firstly consider the case of modes which stay 
under the Hubble radius before the bounce. For that 
modes the initial vacuum state is given by 



-ikrj 



/2k 



(52) 



This can be obtained as a limit of the mode function ((39|) 
for —krj ^ 1. To be specific, let us consider the model 
with a free scalar field and tt^ = O.lZpi. 
At the Hubble radius we have 



Now we are going to perform similar analysis for the 
model with a multi-fluid potential. We obtain formula 



1 



a\H\ 



(53) 



'off 



where 



Pc-n-'piil + wYiSw ~ 1)(1 + 6nGpcil + w)^i^)" x 



which is shown in Fig. [3] We see that for initial time, let 
us say t = — lOOO^pi, all modes with k > 0.003 are well 



4 A72 

3 (l + w)(3u;- 1) 



,„2(1 



3w) 



3 (1 



w) 



described by the function ([52)) . This solutions however 
do not hold during the phase of bounce. Close to the 
(45)bouncc one can approximate 



2 2 + 3u; 



In the limit It I 



3 1- 
00 we obtain 
o 6w - 



2 1 



'off 



(1 



(46) 



(47) 



where we changed time to conformal. Advanced and nor- 
malised solution of the equation in the considered 
limit is 

uik,v) = ^/^JZeMM + i)Hl'^J{-^k) (48) 



'cff ' 



mlsit = 0) = 



1 



(54)1/3 



2/3 _ 



/vg. 



(54) 



For the considered conditions we obtain fco ~ 1.12. In 
this approximation the solutions during the bounce phase 
are 

A, 



Uhc 



B, 



2^1 



'2n 



(55) 



where J7 = ■^/P + fcJ. Finally, in the post-bounce phase 
we have a superposition of advanced and retarded modes 



Mout 



2k 



2k 



(56) 
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FIG. 3: Evolution of the Hubble wave number kn — a\H\ 
Here tt^ = O.lipi. 



In order to derive formulas for the coefficients ak and 
/3fc, we define the matrices 



Mo 



Ml 



M-j 



Ma 



— 7, / flp-'if2'?- 



inri. 



2 



V2k 
k ikrj^ 



\/20 



2 



/2k 



k ikrj^ 



(59) 
(60) 
(61) 
(62) 



Here the relation — |/3fep = 1 holds, as a conse- 

quence of the normalisation condition. Now we have to 
match solutions from the three considered regions to de- 
terminate coefficients and (3^- In order to do that 
we must specify a time when the matching is performed. 
We choose it in the mirror points — 1_ = where 
reaches its maximal value. Then —t- = = to where 



Then matching conditions can be economically written 
as 



Ma^^MaMi^^Mo ' ^ 



Multiplying these matrices we obtain 



(63) 



With use of equation ([M)) we obtain 

2^ 1 L2 ' 6 ' 2 ' -^J 



,1/6- 



For the considered setup we obtain 770 ~ 0.285. 



(57) 



(58) 



—i cos(277ofc) + sin(277oA:) 
at = X 



(ik 



X [2imcos{2riQ9,) + {k^ + fl^) sin(2ryo^^)] (64) 



2kn 



The resulting square of the amplitude for the out state 
modes is 



Uout 



k^ cos[4r/of^] - k'^ sin[2ryof^] {{k + fi)^ sin[2fc(77 - 770) + 2f2?7o] - {k - f2)^ sin[2fc(?7 - 770) - 2^2770] 



r 



(66) 



Based on this result one can calculate power spectrum 
of perturbations. We show this spectrum in Fig. |4l 
Obtained spectrum exhibits sub-horizontal oscillations. 
This effect can be intuitively understood when anal- 
ogy with Schrodinger equation is employed. Namely 
the mode equations are equivalent to a one dimensional 
Schrodinger equation with potential V — —m'^g. Here 
spatial variable is replaced by the conformal time rj. In 
the employed approximation potential is square well of 
width 2ryo and depth 77icg(i — 0). Therefore the evolution 
of the given mode can be seen as transition of a parti- 
cle over the potential well. Amplifications of the ampli- 
tude of transmission correspond to resonances between 



the width of the potential and the phase shift. 

It can be shown that the obtained coefficients {ak,(3k) 
are in fact the Bogolyubov coefficients ak = B+ and 
/3fc = _B_. Therefore the number of produced gravitons 
is given by 



nk = lAf = 



fc^sin^ (2r,ov/^ 



fc2(fc2 







fc2 



(67) 



We show this dependence in Fig. [5] Now it is straight- 
forward to calculate the parameter Jlgw. We show this 
function in Fig. [SI We compare it with the obtained 
low energy approximation f2gw cx; Obtained values 
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FIG. 4: Oscillating tensor power spectrum of the sub- 
horizontal modes at f = 50 Straight line represents spec- 
trum of reference Vt oc . 




FIG. 5: Occupation number of the gravitons in the post- 
bounce state. 




FIG. 6: Parameter f2g„ (bottom curve). Straight line repre- 
sents low energy approximation f2g„ (xv~^. 



of rigw are many orders of magnitude below the present 
threshold for detection. The results obtained were per- 
formed for some simplified model and for the fixed value 

2/3 

of TT^. In particular in the low energy limit fig^ ~ tt^' , 
therefore effect of varying tt^ is considerable. However, 
we do not expect significant changes due to the approxi- 



mations performed. This statement will be confirmed by 
the numerical simulations in Section IVl 



D. Super-horizontal solutions 

In the previous subsection we have shown solutions of 
the mode equation (|13p in the sub- horizontal limit. Now 
we are going to study the super-horizontal A; — )■ limit. 

We introduce new variable in the form 



/ 



(68) 



and change the conformal time to the coordinate one dt — 
adrj. Then equation p3p can be rewritten in the form 



+ n^{k,t)f = 0. 



dt^ 

Here the parameter il'^{k,t) is defined as follows 

Tq 3 / d 1 



n\k,t) 



+ e- 



(69) 



(70) 



We have introduced here parameter the e to trace effects 
of the holonomy corrections in the later equations. In 
the classical limit we should take e = while in presence 
of the holonomy corrections e = 1. Taking k = and 
introducing new complex variable z G C, 



1 



1 



- ^ + i-V6^^c{l + w)t 
we can rewrite equation (|69p in the following form 



^ + Q(.)/ = o 



where 

Qiz[ 

The coefficients are 

ao = 
ai = 

0-2 = 



-I- a\z -f ao 

Z2(z- 1)2 ■ 

9(1 + 2u;) - 4e(l + 3w) 
36(l + u;)2 



(l + u;)2' 
w 



(71) 

(72) 
(73) 

(74) 
(75) 
(76) 



and it will be useful later to remember that ai + a2 = 0. 
Now introducing the new variable 



/(z) = z^(z-l)^.g(z) 



with 



L 
K 



c 
2 

a + b 



(77) 

(78) 
(79) 



one can rewrite equation (|72p as a hypergeometric equa- 
tion 



'Reh 



b+l)z]^-abg = 0. (80) 
dz 



Solution of this equation is given by the hypergeometric 
functions 



g{z) = C2Fi{a,b, c; z). 



(81) 



Furthermore, we have a system of equations for the co- 
efficients 



a2 — afe - 



ao + L{L-l) = 0, (82) 
ai+ab-2KL-2L{L-l) = 0, (83) 
2KL + L{L-1) + K{K -1) = 0. (84) 



One can find that, since ai + a2 ~ 0, we have either 
K ^ L or K ^ 1 - L, where L = i (l ± VI - 4ao) . For 
K = L wc find 



1 



a = - (2c - 1± x/r 
2 ^ 



4Q;ij 

6 = 2c - a - 1, 

c = 2L = 1 ± VI - 4ao, 

and while ii' = 1 — L we have 
1 



a = - (1± Vl + 4ai) , 



1 — a, 

2L = 1 ± vr 



4ao 



(85) 

(86) 
(87) 

(88) 

(89) 
(90) 



As an exemplary solution we consider w — 1 case both 
classically (e — 0) and with holonomy corrections to the 



mode equation (e 
we have 



where 



c±(e=l) = l± 



1). Then since ai = —a^ 



— and c+ (e = 0) = 1 ± 

6 ±V 7 2 



-1/4 
(91) 

(92) 



In Fig. [7] we show solutions for the real components of 
h variable. 

In Fig. [8] we show solutions for the imaginary compo- 
nents of the h variable. 

In Fig. [niwe show solutions for the absolute value of 
the h variable. 

As it can be seen, solutions with and without the quan- 
tum holonomy corrections are qualitatively similar. An- 
other observations is that for the times t ^ 1 evolution 
takes a logarithmic form. This result can be understood 
considering equation (|13p in the classical limit and taking 
k —^ 0. Then one can find approximate solution in the 
form 



Ak + Bi 



_drf_ 



(93) 




FIG. 7: Real components of h variable. Dashed (blue) line 
represents solution with c+(e = 1) while straight (red) line 
represents solution with c+(e — 0). 

Jmh 




<1 



FIG. 8: Imaginary components of h variable. Dashed (blue) 
line represents solution with c+(e = 1) while straight (red) 
line represents solution with c+(e = 0). 




FIG. 9: Absolute values of h variable. Dashed (blue) line 
represents solution with c+(e = 1) while straight (red) line 
represents solution with c+(e — 0). 



where Aj. and Bk are some constants. For the models 
considered in the present paper we have a cx (±7y)^/ (i+3io) 
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and ±r] oc {±t)T^(^ 



panding phase and 
fore for the considered w 



. Here we have + sign for the ex- 
sign for the contracting one. There- 
1 case we find 



/i ~ ife + Bk \n{±t) for \t\ > 1 (94) 
in agreement with the solutions found in this subsection. 

V. NUMERICAL INVESTIGATIONS 

As it was shown in the previous section, analytic so- 
lutions of the mode equation are available only in some 
limits. Namely for both t'^ ^ oo and fc ^ 0. Also for 
fc — !■ oo an approximate solution was found. It is however 
not sufficient to describe whole spectrum of the gravita- 
tional waves produced on the bounce phase since the in- 
teresting intermediate regimes are unexplored. Therefore 
numerical analysis is required. 

In the numerical computations we are going to solve 
the autonomous system of equations 



du 
dr] 
d-iTu 
dr] 
dt 
drj 



(95) 
(96) 
(97) 



where a{t) and 'Ti^g(i) are defined for particular back- 
ground dynamics. In the considered models with free 
scalar field and multi-fluid potential these functions are 
given by analytical expressions. Since canonical variables 
M, 7r„ g C we decompose 



U = Ul + iU2: 



(98) 
(99) 



Now it is crucial to define proper initial conditions for 
(mi, M2, 7r„i, 7r„2) for some time r^g. It is always unam- 
biguous how to choose a proper vacuum defined on the 
cosmological background. However on the sub-horizontal 
scales, when Minkowski space approximation holds, we 
can set 



and 



7I"«i(?7o) 



1 



cos(fc77o), 
sin(fc77o) 



-\l-sm{kT]o) 



- cos(fc77o) 



(100) 
(101) 

(102) 
(103) 



at some time ryo- Here we set initial values like in the 
model of sub-horizontal modes studied in the previous 



section. Therefore analysis is correct for the modes with 
with k > 0.003. 

In Fig. [TO] we plot evolution of the A: = 0.1 mode 
during the bounce phase. We compare here the evolu- 




FIG. 10: Evolution of the modes with k = 0.1. Dotted (blue) 
curve represent solution of mode equations with holonomy 
corrections. Straight (red) curve represent solution of mode 
equations without holonomy effects. 

tion of modes with and without holonomy corrections to 
the modes equation. We see that close to the turning 
point the effects of the holonomy corrections become sig- 
nificant. However the further oscillating evolution does 
not change qualitatively. The difference is some suppres- 
sion of the amplitude of perturbations due to the quan- 
tum corrections. This feature can be also seen in Fig. 
111! where classical and quantum corrected tensor power 
spectra are shown. We find, comparing with the clas- 



with holonomy corrections 




t=50[/pi] 



FIG. 11: Comparison between sub-horizontal spectra ob- 
tained with and without holonomy corrections to mode equa- 
tion. Dashed (blue) line represents rescaled initial vacuum 
power spectrum. 

sical case, that quantum holonomy effects amplify low 
energy modes. Therefore tensor power spectra increases 
by about one order of magnitude. For the high energies 
classical spectra starts to dominate slightly. It is also 
worth to notice that oscillations do not overlap. 
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To impose initial conditions on the super-horizontal 
scales one can use instantaneous vacuum. This is how- 
ever possible only for values of k fulfilling > 0. As we 
have found earlier this condition is fulfilled for all k in the 
model with the free scalar field. Therefore initial instan- 
taneous vacuum state can be defined on all length scales. 
It can be shown that Hamiltonian of perturbations at 
time 770 is minimised for 




part becomes nearly flat (depending on the model of in- 
flation). It is possible that the oscillating features will 
also survive giving the footprints of the bouncing phase. 
However the further analysis has to be performed to ap- 
prove these speculations. In particular inflationary power 
spectrum with the obtained bouncing initial conditions 
must be calculated. Then it will be possible to compute 
the B-type polarization spectra of CMB. Therefore a way 
to relate the quantum cosmological effects with the low 
energy physics becomes potentially available. 



VI. SUMMARY 



where cuk = \/k^ + m^g. 

In Fig. [T3 we show tensor power spectrum at post- 
bounce stage (t = 50 Ipi) with imposed instantaneous 
vacuum initial conditions at t = —1000 Ipi. The charac- 




FIG. 12: Whole post-bounce tensor spectra in presence of the 
holonomy corrections. Green points comes from the numerical 
simulations. Black line is the analytical spectrum from the 
model given by Eq. 1661 Dashed (red and blue) lines represents 
UV and IR behaviours, in both cases Vt oc . 

teristic feature of the spectrum are oscillations. Moreover 
both UV and IR behaviours are in the form Vt oc k^ . 
We see that analytical model given by Eq. [66] fairly good 
overlap with the numerical results. Especially structure 
of oscillations is exactly recovered. Also the asymptotic 
behaviours are consistent. The evident discrepancy is the 
difference in the total amplitude. In fact this difference 
can be suitably adjusted varying parameters of the model 
?7o and fco. Then low energy behaviour can be exactly re- 
covered. However it introduces additional phase shift and 
structures of oscillation no longer overlap. It is also im- 
portant to note that effect of the imposed instantaneous 
vacuum initial conditions is negligible in the range stud- 
ied. Therefore Minkowski vacuum approximation is still 
valid. 

The obtained power spectrum can be now applied as 
an initial condition for the inflationary modes. We ex- 
pect that super-horizontal part of the spectrum does not 
change during the inflationary phase. However the UV 



In this paper we have considered influence of Loop 
Quantum Gravity effects on the gravitational waves prop- 
agation in the flat FRW cosmological background. The 
considerations presented based on the semi-classical ap- 
proach where quantum effects are introduced by correc- 
tions to the classical equations of motion. This approach 
was successfully applied to the homogeneous models. In 
this case good agreement between results of fully quan- 
tum and semi-classical analyses was found. Here we have 
applied semi-classical approach to the inhomogeneous 
model where inhomogeneity is treated perturbatively. 
Therefore perturbations had no influence on background. 
In general both classical and quantum back-reaction ef- 
fects can be important close to the phase of bounce. Here 
we assumed that they can be neglected. We have also not 
considered effects of the quantum fluctuations of back- 
ground on the inhomogeneities. Quantum effects were 
introduced by the so called holonomy corrections. In the 
homogeneous models these corrections lead to absence 
of the initial singularity and emergence of the bounce 
phase. Effects of the other known type of LQG correc- 
tions, the inverse volume ones, were studied earlier in nu- 
merous papers. Here we considered self-consistent model 
where holonomy effects influence both background and 
perturbations (gravitational waves). In the earlier stud- 
ies effects on background and perturbations were studied 
independently. In particular in Ref . [l^ a model of gravi- 
tational waves production during the holonomy-inducted 
bounce phase was investigated. In Ref. 15l,117|] effects of 
holonomy corrections on the gravitational waves in infla- 
tionary phase were studied. However quantum effects on 
the background dynamics were neglected there. More- 
over quantum-corrected source term was not taken into 
account in those studies. Linear part of this term van- 
ish in the classical limit. However, its contribute while 
holonomy corrections are present. Therefore source term 
has to be taken into account in the full treatment. In the 
present paper we have included effects of this term. 

We have considered models with both free scalar fleld 
and self-interacting fleld with multi-fluid potential. In 
both cases scalar fleld is a monotonic function and can 
be treated as a internal time variable. 

We have shown that in the model with the free fleld, 
effective mass term m^g for gravitational waves is a non- 
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negative function. This is not the case for the models 
with multi-fluid potential. We have found solutions of 
the mode function in the pre-bounce phase and deter- 
mined the power spectra of the obtained perturbations. 
Then we have considered sub-horizontal solutions dur- 
ing the bounce phase. We matched the solutions from 
pre-bounce, bounce and post-bounce phases. Based on 
this we have found power spectrum of gravitational waves 
and determined Bogolyubov coefhcients. Then number 
of produced gravitons rtk and the parameter fig^ were 
calculated. We have found that fig^ reaches 10~^^ which 
is far below any observational bounds. These results were 
obtained for fixed parameter tt^ = 0.1 Zpi. 

Based on analytical considerations we have found that 
power spectrum exhibits oscillations on sub-horizontal 
scales. An intuitive explanation of this effect was given. 
We have also solved the model analytically in the super- 
horizontal limit. These results indicate that quantum 
corrections do not introduce qualitative difference in the 
power spectrum on these scales. Therefore the obtained 
lack of power on the large scales is a feature of the bounc- 
ing evolution and not of the quantum corrections to the 
mode equation. 

Subsequently we have investigated the model numer- 
ically. We have approved presence of the oscillations 
emerged from the simplified analytical considerations. 
Both numerical and analytical results were compared. 
We have found good qualitative and quantitative agree- 
ment. We have also approved earlier observation that 
quantum corrections does not introduce qualitative dif- 
ference in the power spectrum. The only differences ob- 
served were in total amplitude and phase of oscillations. 

Imposing initial instantaneous vacuum state we have 
also studied the low energy part of the power spectrum. 
Therefore we have finally found the full shape of the ten- 
sor power spectrum. This spectrum can be used to study 
further phenomenological consequences. In particular, it 
can be applied as an initial condition for the inflationary 
modes. Then we expect that the sub-horizontal part of 
the spectrum becomes flat while super-horizontal form 
survive. It is also possible that sub-horizontal oscilla- 
tions survive as features of the dominant nearly flat in- 
flationary spectrum. Therefore two observational effects 
of the bouncing phase can be distinguished: oscillations 
and lack of power on the super-horizontal scales. These 
effects can potentially be tested with the future CMB 
missions like Planck H or proposed CMBPol Es- 
pecially promising are observations of the CMB polar- 
ization. Here bounce can lead to the low multipoles sup- 
pression in the B-type spectrum. At present projects 

27[ are aiming 



like Clover |25l, QUaD [26| or QUIET 
to detect this spectrum and first results are expected in 
the near future. Therefore it is the next step to derive 
quantitative predictions of the CMB features from the 
presented model. 
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APPENDIX A: LOOP QUANTUM GRAVITY 
WITH GRAVITATIONAL WAVES 



Loop Quantum Gravity (LQG) describes the gravita- 
tional field as SU (2) non-Abelian gauge field using back- 
ground independent methods. The canonical fields are 
so called Ashtekar variables {A = A^ndx'^.E = Efr^da) 
which take value in bu(2) and su(2)* algebras respectively 
and they fulfil the Poisson bracket 



(Al) 



where k = SnG and 7 is the Barbero-Immirzi parame- 
ter. These variables are analogues of the vector potential 
and the electric field in electrodynamics. The Ashtekar 
variables are related with triad representation. In LQG 
gauge fields describe only spatial part S when time is 
treated separately. 

In cosmological applications we perturb basic variables 
around a background 



At = Al- 



SAl. 



(A2) 
(A3) 



For the spatially fiat FRW background components have 
the following form 



e: = ps:, 



(A4) 
(A5) 



where p = a? and k = p/2p. Perturbations can be split 
for the scalar, vector and tensor parts. For the purpose of 
this paper we consider here only the gravitational waves 
(tensor part). Tensor perturbations of the flat FRW met- 
ric are introduced as follows 

500 = -N^ + qabN'^N'' = -N = -a'', 

gab = Qab ^ a^[Sab + hab], 

with the conditions = dah'^ = and \hab\ ^ 1. In the 
TT gauge h\ = — /ij = ft.® and h2 — ~ h^. 

Now we are going to perturb the Hamiltonian of the 
theory. The full Hamiltonian is composed of the gravi- 
tational and Hq and matter parts. Hamiltonian Hq 
takes the form of a liner combination of the constraints 



Hr 



d^^iN'O^ + N^Ca + NS). 



Spatial diffeomorphisms constraint: 

Ca = E^F:, - (1 - ^^)KG. 
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Gauss constraint: 
Scalar constraint: 



5* = 



'feF^,-2(l + 7^)iffX] 



where F = dA + \ [A, A\. However thanks to the quan- 
tum gravity effect this Hamiltonian undergoes modifi- 
cations. These modifications can be introduced on the 
phenomenological level by the replacement 



k 



sin nfijk 
nflj 



(A6) 



in the classical expressions. Here 



where A — l^i-K^l 



This kind of corrections we call holonomy ones. Factor 
n can be fixed from requirement of the anomaly cancel- 
lation [3j 113 ■ Effective second order Hamiltonian with 
holonomy corrections takes the form 



rrphen -1 



d^xN 

s 

sin p.'jk 



sin iijk 



2p3/2 Y 

TP (^) 



+ ^iS.,6^'6^fd.E^djE, 



where for tensor modes 
5Et ^ 

Based on the Hamilton equations 



1 


K + 


( sin2^7fc^ ^^ 


2 




\ 2a27 ) ^ 



H. 



H: 



SEf = [SEf, 

5ki = [ski ^ 

we obtain equation 
where 



phcn 
G 



phcn 



H„ 
H„ 

167rGH: 



Qa 



Tq 



.Ai dp. 



sin ^jk 

Pi 



(A7) 

(A8) 
(A9) 

(AlO) 
(All) 

(A12) 
(A13) 



1 dH^ 
Wq dp 



— ^ cos2/i7A;- 



SH, 



6iSE^ 



4) 



are quantum holonomy corrections. 

We consider homogeneous scalar field with the Hamil- 
tonian 



where up to the second order 

1 



^\dctE\V{^) 



(A15) 



Vdet^ 



p2 



VdeiTE 



1 



p2 



-^5l5E^5i5E] 



(A16) 



±,6:8E<;5l6Et 



(A17) 

^ SiSE'' = the above ex- 



pansion simplifies. Then 



i/,„ = / ^ I 1 ^ _ y^^A siSE^SiSE'^+0{E^) 



Now we can derive variation 



SH„ 



1 N 



5{5Ef) 2Vp\2p3 
and derivative 



V{<i^) SlS'JEt (A18) 



(A19) 



One can now easily find that in the classical limit, 
when we set cos{2 fljk) = 1 in expression (|A14p . the 
source term vanish. This is due to the opposite signs of 
the bracketed expression in equations (jAlSp and (jA19p . 
When quantum holonomy corrections are present we have 



cos(2/27fc) = 1-2—, 

Pc 



(A20) 



which can be found from background equations of mo- 
tion. Therefore the form of the quantum corrections sim- 
plifies to 



Tc 



HI 



SttGpp^ 
3 Pc ' 

^ Pc 



(A21) 
(A22) 



where we have chosen N — ^Jp and adopted the expres- 
sion (IMI). 
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